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Gutzwiller variational method is applied to an effective two-dimensional Hubbard model to ex- 
amine the recently proposed gossamer superconductor by Laughlin [[i]]. The ground state at half 
filled electron density is a gossamer superconductor for smaller intra-site Coulomb repulsion U and 
a Mott insulator for larger U. The gossamer superconducting state is similar to the resonant valence 
bond superconducting state, except that the chemical potential is approximately pinned at the mid 
of the two Hubbard bands away from the half filled. 



Theories for high temperature superconductivity con- 
tinue to attract much interest in condensed matter 
physics. Soon after its discovery, Anderson proposed 
the idea of resonant valence bond (RVB) state for the 
observed unusual properties in high T c superconducting 
Cu-oxides[||. In the RVB picture, each lattice site is ei- 
ther unoccupied or singly occupied by a spin-up or down 
electron. The spins are coupled antiferromagnetically 
without long range order. The charge carriers move in 
the spin liquid background and condense to a supercon- 
ducting state(|, |j. The RVB states are oftenly studied 
using 2-dimensional Hubbard or t- J models ^ |i| . In this 
scenario, the undoped cuprate with density one electron 
per site, or half filled, is a Mott insulator, and the su- 
perconductor is viewed as a doped Mott insulator when 
additional holes or electrons are introduced. A Mott in- 
sulator is a special type of insulator casued by electron 
interaction. It has been established that the ground state 
of many models in 1-dimensional chain or in ladders at 
half filled are Mott insulator ||. In two or higher dimen- 
sions, a Mott insulator has a strong tendency toward an- 
tiferromagnetic or other types of ordering states breaking 
translational symmetry. 

Very recently, Laughlin has proposed an interesting 
new notion, the gossamer superconductor, for high tem- 
perature superconducting Cu-oxidesQ. In a gossamer 
superconductor, the superfluid density is very thin, in 
contrast to the conventional superconductor. Laughlin 
has proposed an explicit microscopic wavefunction for the 
gossamer superconductor, which has the following form, 
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where | ^bcs) is a BCS superconducting state, and 
II Q = JX(1 — ocrii^nii) is a projection operator to par- 
tially project out doubly occupied electron states on each 
lattice site i. ni a = c^Cfo is the electron number op- 
erator of spin a =1, j at site i, and a is a parame- 
ter between and 1. If a is close to 1, the projection 



operator strongly suppresses the superfluid density. In 
the gossamer superconducting state, the probability to 
have both spin-up and spin-down electrons occupying the 
same lattice site is largely reduced but remains to be fi- 
nite. Because of the partial projection, the state is super- 
conducting even at half filling. This is different from the 
RVB theory where the electron doubly occupied states 
are completely projected out, hence the half filled RVB 
state is a Mott insulator and the superconductivity oc- 
curs only away from the half filled. As it has been shown 
by Laughlinjj, the gossamer superconducting state in 
Eq. (1) is an exact ground state of a model Hamilto- 

nian given by H L = Y.ka E k}\}k^ where E ka ^ °' 
bt = n a bz n„ , and is the quasiparticle annihila- 
tion operator of the BCS state, bjg- | ^bcs) = 0. n^ 1 
is the inverse of II Q . Laughlin has also argued that the 
gossamer superconducting state is related to the large 
on-site Coulomb repulsion. 

It will be interesting to examine the possible gossamer 
superconducting state in a more realistic model. As it is 
generally believed that the large Coulomb repulsion may 
lead to a Mott insulator at half filling, it will also be 
interesting to examine the possibility of the phase tran- 
sition from a Mott insulator to a gossamer superconduc- 
tor as the electron interaction strength decreases. Since 
a doped Mott insulator can be a RVB superconducting 
state, it is natural to ask the question of the similarities 
and the differences between the gossamer and the RVB 
superconductors. 

In this paper we intend to examine these questions by 
studying an effective Hubbard model given in Eq. (2) 
below using Gutzwiller's variational method R|. In the 
Gutzwiller's approach, the on-site Coulomb repulsion is 
treated exactly, while the kinetic energy is studied vari- 
ationally, so that it is suitable to examine some issues 
in strongly correlated systems. That method was used 
by Brinkman and RiceQ to study the phase transition 
between an insulator and a metallic state described by 
a partially projected Fermi liquid state. The variational 
method applied to the effective Hubbard model in two- 
dimension demonstrates a phase transition from a gos- 
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samer superconductor for smaller intra-site Coulomb re- 
pulsion U to a Mott insulator for larger U at half fill- 
ing. The gossamer superconductor is shown similar to the 
RVB superconducting state of the doped Mott insulator. 
However, its chemical potential is found to be approx- 
imately pinned at the mid of the lower and the higher 
Hubbard bands, different from the RVB state where the 
chemical potential is shifted to the lower Hubbard band 
upon doping. 
We study an effective Hubbard Hamiltonian, 



H 



+ h.c.)+J2 JijSi-Sp) 

(ij) 



In this Hamiltonian, we have introduced an antiferro- 
magnetic spin-spin coupling term (Jy > 0) to account 
for the virtual electron double occupancy effect. In the 
large U limit, rs 4t?-/£7. This model may be viewed 
as an effective Hamiltonian of the Hubbard model. The 
inclusion of the antiferromagnetic spin coupling appears 
consistent with the weak coupling renormalization group 
analyses [||, and is appropriate in the variational ap- 
proach studied here. Although the precise values of Jy 
are to be determined, that does not alter the qualitative 
physics we will discuss in this paper. In the limit U — > oo, 
the model is reduced to the t-J model. 

We consider | in Eq.(l) as a variational trial wave- 
function to examine the superconductor- insulator tran- 
sition at half filling, and to compare the gossamer super- 
conducting state with the RVB state away from the half 
filled. In our theory, tig, v^, and a are variational param- 
eters. In the limiting case upj^ — 0, | ^bcs) reduces to 
the Fermi liquid state, a = corresponds to the uncor- 
rected state (U — 0), and a = 1 corresponds to the limit 
of no doubly occupied state. Obviously a = 1 if U — > oo. 

The variational energy E = (H) is given by, 



E = Ud+(H t ) + (Hj) 



(3) 



where d = (n^n^) is the electron double occupation 
number, and < d < 1. (A) is the expectation value (per 
site) of operator A in the state | The first term is 

the intra-site Coulomb interaction energy, while the sec- 
ond and the third terms are the average kinetic and spin- 
spin correlation energies, respectively. Note that in the 
Gutzwillcr approach and at half filling, d is a measure of 
the mobile carrier density n c and proportional to n c /m* 
measured in the a.c. conductivity with m* the effective 
mass. At d = 0, i.e. a = 1, we have (H t ) = 0, and 
E = (Hj) at the half filling. This state describes a Mott 
insulator. The case with d>0or0<a<I describes 
Laughlin's gossamer superconducting state. 

The variational parameter d is determined by the con- 
dition dE{d)/dd = 0, or 



U + d(H t ) /dd + d(Hj) /dd = 



(4) 



At half filling, we expect a transition from the Mott in- 
sulator at larger U to the gossamer supercondutor at 



smaller U as U decreases passing through a critical point 
U c . To study this phase transition, we follow Brinkman 
and Rice || and compare the energies of the two states 
with the difference that here we consider the projected 
BCS state while Brinkman and Rice considered the pro- 
jected Fermi liquid state. The transition point U c is given 
by U c = {-d(H t ) /dd-d(Hj) /dd) u=a . For U > U c , there 
is no solution of Eq. (4) for physical values of d, indicat- 
ing that d = 0. 

We use the Gutzwiller approximation B to carry out 
the variation and to estimate U c . In the Gutzwiller ap- 
proximation, the effect of the partial projection operator 
on the doubly occupied sites is taken into account by a 
classical statistical weighting factor which multiplies the 
quantum coherent result calculated for the unprojected 
state | ^bcs)- A clear description of the method has 
been given by Vollhardt |T(|. The method was used to 
study the 2-dimensional t-J model where the pro- 
jection operator corresponds to the case of a = 1 (the 
complete projection). In the present model, the hopping 
and the spin-spin correlation energies in the state \ s S>l) 
are related to those in the unprojected state | ^ bcs) by 
the corresponding renormalized constants gt and g s : 



(cLoj 



(Si ■ Sj) = g s (S l ■ Sj} 



(5) 



where (A)q is the expectation value of operator A in 
the state | ^bcs)- The renormalization facotrs gt and 
g s axe determined by the ratios of the probabilities of 
the corresponding physical processes in the states | ^l) 
and | ^>bcs)- By counting the probabilities we ob- 
tain these renormalization constants for the partially pro- 
jected state (n: electron density), 
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(n - 2d)(Vd + y/l-n + d) 2 
(1 - n/2)n 

r (n-2d) a 
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(6) 



The value of gt is the same as that previously obtained 
for the projected Fermi liquid state [l0|]. At d = 0, we 
have g t = 26/(1 + 5) and g s = 4/(1 + S) 2 , with 6=1- n, 
recovering the results in Ref. |Uj. At half filling, n = 1, 
we have g t = 8d(l — 2d) and g s = 4(1 — 2d) 2 . Using the 
Gutzwillcr approximation, The variational condition for 
d in Eq. (4) becomes 



U- 



dgt 
dd 



H, 



to 



dg s 
dd 



H 



JO 



= 



(7) 



The transition point between the Mott insulator and 
the gossamer superconductor is thus found to be, U c — 
-8{H t )o + 16(Hj) Q . Note that (H t ) < 0, so that U c is 
generally positive if the kinetic energy term in the uncor- 
related state dominates. . In the insulating phase, only 
the spin-spin interaction is non-zero. The problem be- 
comes identical to that in the RVB theory at half filling, 
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and there is a redundancy in the fermion representation 
of the state due to a local SU(2) symmetry of the spin 
Hamiltonian [O, The redundancy is removed in the 
gossamer superconducting state for the kinetic energy 
term breaks the SU(2) symmetry, similar to the effect 
of doping in the t-J model. At d << f, the symme- 
try of the gossamer superconductivity is the same as the 
symmetry of the RVB state. Within the Gutzwiller ap- 
proximation, the pairing order parameter in the gossamer 
superconductor is related to the uncorrelated state by a 
renormalized factor g t , 



(8) 



Near the transition point, g t = 8d « 1, indicating the 
smallness of the superfluid density, a quantitative mea- 
sure of the gossamer superconductivity. It is interesting 
to note that the pairing order parameter in the RVB state 
has also the form of Eq. (8) with gt = 26 for 6 « 1. 
This comparison indicates that a gossamer superconduc- 
tor with double occupation d at half filling is similar to 
the RVB superconductor at doping S with the correspon- 
dance of 6 = id. 

In what follows we take an example and consider the 
effective Hamiltonian in a 2-dimensional square lattice 
with only the nearest neighbour hopping iy = t and the 
nearest neighbor spin coupling Jy = J and consider the 
case n < 1. For any given value of d, the Coulomb inter- 
action term in the present theory contributes a constant 
Ud to the variational energy, and the variational proce- 
dure for other parameters (u^ and v-A is almost the same 
as that in study of the t-J model carried out in Ref. |ll[] 
except that the renormalization constants gt and g s here 
depend also on the double occupation d. 

We introduce two correlation functions (r = x,y), 

A r = E^CifCi+T-^O, Xt = E,x (4t C *+t i( t)o- TIlC 

variational solution is then given by the coupled gap 
equations, 

k 

where = y] r A T cosfc r , = z% ~ Z) T Xt cosk r . 



In the above equations, = y| | 2 +Xp e "g — 

[— 2g t t(cos k x -\- coskyj /2]/(3(7 S t//4), and jl is related 
to the chemical potential n by 



(i = fi+ —{H t )o + -^-{Hj) 
On On 



(9) 



In Eq. (9), the second and the third terms originate 
from the n dependences of g t and g s in the variational 
procedure pT| , which will be important in calculation of 
the chemical potential of the state. These gap equations 



must be solved simultaneously with Eq. (7) for d and the 
electron number equation given by S = Y^kXg/^g- 

We first discuss the half filled case. The ground state 
of the insulating phase (d — 0) is the same as that of the 
Heisenberg model. In the metallic phase (0 < d « 1), 
the kinetic energy breaks the local SU(2) symmetry and 
favors the d-wave superconducting state with A^,- = — A y . 
The symmetry is the same as the symmetry studied in 
the t-J model slighlt away from the half filled Jll], [ll| 
|T^| . At the superconductor- insulator transition point, 
we have (H t } = -2^2tC, and (Hj) = -(3/4)JC 2 , 
with C — \ \J cos2 k x + cos 2 k y = 0.479. We estimate 
from these values that U c = 10. 8t — 2.75 J. 

We now discuss the slightly less than half filled case. 
We expect that the variational parameter d is a smooth 
function of the electron density around the half filled. 
The gossamer superconducting state essentially remains 
unchanged in the regime 5 << d, and the superconduct- 
ing order parameter is mainly controlled by d, weakly 
depending on 5 as we can see from the expression for 
gt- The chemical potential [i can be calculated by us- 
ing Eq. (9). In the limit 5 — ^ + , /i — > 0, and we have 
\i -► -4(1 - 4d)(H t )a + 8(1 - 2d)(Hj) = U/2. In the 
last step of the above calculations, we have used varia- 
tional equation (7) to relate the kinetic and spin coupling 
energies to the Coulomb energy. Since \i = U/2 at the 
half filled by electron-hole symmetry of the model, we 
conclude that in the gossamer superconducting state the 
chemical potential is continuous at the half filled, and is 
pinned at the mid of the lower and the higher Hubbard 
bands. This result is reasonable because the gossamer su- 
perconducting state is a metallic state and the chemical 
potential is expected to be continuous 0. This feature 
is in contrast to the RVB state discussed below. 

Ii U > U c with U c defined as the critical U at half 
filling, the state changes dramatically from an insulator 
to a RVB superconducting state as the electron density 
varies away from the half filled. At U » U c , d changes 
very little from zero the physics is essentially the 
same as that given by the t-J model. While the RVB 
state is similar to the gossamer superconducting state 
in the sense that they have the same pairing symme- 
tries and small pairing order parameters, the chemical 
potential in the RVB state is very different from that in 
the gossamer superconductor. To see this explicitly, we 
consider the limit 5 — > + , so that we have fx — ► 0, and 
dgt/dn — — 2, dgj/dn = 8. From these values, we obtain 
H -> -2(H t )o + 8(Hj) ~ 2.7t - 1.38J, which is about 
U c /4: for J « t, and is much smaller than the chemi- 
cal potential U/2 at the exact half filled. We conclude 
that the chemical potential in the RVB state is discon- 
tinuous at the half filled, and it is shifted from U/2 at 
the half filled to the lower Hubbard band away from the 
half filled. The difference in the chemical potentials in 
the gossamer and the RVB states can in principle be dis- 
tinguished in spectroscopic experiments, although other 
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symmetry broken states not included in the variational 
theory and the inhomogeneity will complicate the analy- 
ses.. 

The actual ground state of the Hubbard model in a 
square lattice and at half filled is an antiferromagnet, 
which we have not included in our variational wavefunc- 
tion. Nevertheless, the phase transition and the sim- 
ilarities and the differences between the gossamer and 
the RVB superconductors should be relevant to the sys- 
tems away from half filled. It is interesting to note that 
the half filled superconducting state may be stabilized 
against the antiferromagnetism in the presence of an ex- 
plicit pair hopping term in the Hubbard model as stud- 
ied by Assaad et al. Jl6| . The gossamer superconductor 
should also be relevant to the frustrated magnetic sys- 
tems where the anitferromagnetism is suppressed. For 
instance, there have been numerical studies of the an- 
tiferromagnetic Heisenberg spin model with nearest and 
next nearest neighbor couplings to show evidences for a 
spin liquid ground state [[l7|. There have been numeri- 
cal studies of the Hubbard model at half filled with the 
nearest and next nearest hoppings to demonstrate a para- 
magnetic insulating phase with a transition to a metallic 
phase at zero temperature as U decreases |lj, It will 
be interesting to examine the possible gossamer super- 
conductivity in that metallic phase. 

In summary, we have used Gutzwiller variational 
method to study an effective Hubbard model. The cal- 
culation based on the Gutzwiller approximation supports 
Laughlin's recent proposal of gossamer superconductor at 
relatively smaller infra-site electron Coulomb repulsion 
U, and predicts a phase transition from the gossamer su- 
perconductor to the Mott insulator as U increases at den- 
sity 1 electron per site. The gossamer superconductor is 
similar to the RVB superconducting state with the major 
difference on the positions of their chemical potentials. 
The Gutzwiller approximation we used in this paper has 
been previously tested against variational Monte Carlo 
method [^(j) with quite good agreement [|llj . We believe 
that the qualitative conclusions obtained here should be 
reliable, while the quantitative values should be taken 
with care and may be refined using variational Monte 
Carlo or other numerical calculations. . 
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